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Abstract 

In the framework of the /^-factorization QCD approach we consider the production of 
b quark pairs in pp collisions at the Fermilab Tevatron. We investigate the dependence of 
the b quark, B meson and decay muon differential cross sections on the different forms of 
unintegrated gluon distributions. The analysis also covers the azimuthal correlations between 
the b and b quarks and their decay muons. Our theoretical results agree well with recent data 
taken by the D0 and CDF collaborations at Tevatron. Finally, we present our predictions 
for muon-muon and muon-jet cross sections at the Tevatron and CERN LHC conditions. 



1 Introduction 

Recently D0 and CDF collaborations have reported new experimental data [1-5] on the 
6-flavor production at the Tevatron. These data are found to be about a factor of two or 
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more larger than the predictions of perturbation theory (pQCD) at next-to-leading order 
(NLO) [1-6]. Therefore, it would be certainly reasonable to try a different way. 

At the energies of modern colliders (such as Tevatron and LHC), heavy quark and quarko- 
nium production processes belong to the class of the so called semihard processes [7-10]. In 
these processes, by the definition, the hard scattering scale /i ~ itlq is much larger than the 
QCD parameter Aqcd, but, on the other hand, it is much smaller than the total center-of- 
mass energy: Aqcd <C /x <C y/s. The last condition means that these processes occur in 
the small x region, x ~ itlq/^/s <C 1, where the contributions of the "large logarithms" to 
the evolution of gluon densities are known to become rather important. It is known also 
that in the small x region it becomes necessary to take into account the dependence of the 
subprocess cross sections and gluon structure functions on the gluon transverse momentum 
kx [7-10]. Therefore, the /c T -factorization (or semihard) approach provides a more suitable 
ground for the calculations than the ordinary parton model. 

The /^-factorization approach is based on the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) [11] evolution equation for gluon densities. The resummation of the terms 
a n s ln n (/i 2 /A QCD ), a n s ln"(/j 2 /A QCD ) ln n {l/x) and a n s ln n (l/x) leads to the so called unin- 
tegrated (g^-dependent) gluon distributions $(x,gf,,/i 2 ) which determine the probability to 
find a gluon carrying the longitudinal momentum fraction x and transverse momentum qx 
at the probing scale /i 2 . In contrast with the usual parton model, the unintegrated gluon 
distributions have to be convoluted with off-mass-shell matrix elements corresponding to 
the relevant partonic subprocesses [7-10]. In the off-mass-shell matrix elements, the virtual 
gluon polarization tensor is taken in the form [7-10]: 

If = &&. (1) 
Qt 

The 6-flavor production at the Tevatron in the /c^-factorization approach was considered 
earlier in [8, 12-17]. However, the off-shell matrix elements of hard partonic subprocess 
presented in [8] contain some misprints. Also, the results obtained in the paper [16] do not 
agree with other results [12-15]. 

In our previous paper [17], we analyzed the px distribution of the produced 6-quarks 
(presented in the form of integrated cross sections). Our results agree well with D0 [5] and 
CDF [3] experimental data. Also, we found that our off-shell matrix elements for partonic 
subprocess coincide with the ones presented in Ref. [9]. 

Here, we use the /c^-factorization approach for a more detailed analysis of the experimental 
data [1-5]. We inspect the dependence of the b quark, B meson and decay muon cross sections 
on the different forms of the unintegrated gluon distributions. The analysis also covers the 
azimuthal correlations between the b and b quarks and their decay muons. Special attention 
is paid to the role of the unintegrated gluon distributions which has been applied earlier in 
our previous papers [17-22]. In addition, we present our predictions for muon-muon and 
muon-jet cross sections at the Tevatron and CERN LHC conditions. 

The outline of this paper is as follows. In Section 2 we present the analytical expressions 
for the total and differential cross sections of the inclusive heavy quark production in the 
/c T -factorization approach and describe the unintegrated gluon distributions which we use 
in our calculations. In Section 3 we present the numerical results of our calculations and 
compare them with the D0 [4, 5] and CDF [1-3] data. Finally, in Section 4, we give some 
conclusions. 
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Figure 1: Diagram for the pp — > bbX process. 



2 Theoretical framework 

In this section we present the expressions for the inclusive heavy quark production total 
and differential cross sections in the A/p-factorization approach and describe the unintegrated 
gluon distributions which we use in our calculations. 



2.1 Kinematics 

As indicated in Fig. 1, we denote the 4-momenta of the incoming protons and the outgoing 
heavy quarks as k%, k 2 and pi, p 2 , respectively. The initial gluons have the 4-momenta q\ 
and q 2 . We use the Sudakov decomposition, which has the form 

Pi = a x ki + (3ik 2 +Pit, P2 = a 2 k x + (3 2 k 2 +p 2T , ^ 
qi = xiki + qiT, q2 = x 2 k 2 + q 2Tl 

where pit, p 2 r, Qit an d q%r are the transverse (4-vector) momenta of the corresponding 
particles, and 

2 2 2 2 2 2 2 /o\ 

Pl=?2 = m Qi Ql = QlTi ?2 = ?2T- ( 3 ) 

In the pp cm. frame we can write: 

fci = v^/2 (1, 0, 0, 1), k 2 = (1, 0, 0, -1), (4) 

where we neglect the masses of the protons. The Sudakov variables are expressed as follows: 

miT i \ ™2T , s. 

a x = — =-exp(j/ij, ol 2 = -^exp(y 2 ), 

Vs (5) 

TTli r p Tn 2 T y ' 

Pi = —j=r exp(-j/i), p 2 = —=r exp(-y 2 ), 
V s V s 

where m\ 2T = m 2 Q + p\ 2T , and yi and y 2 are the rapidities of final heavy quarks in the pp 
cm. frame. From the conservation laws, we can easy obtain the following conditions: 

xi = at + a 2 , x 2 = Pi + (3 2 , qvr + <l2T = Pit + P2T- (6) 
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2.2 Inclusive heavy quark production cross section 



Here we recall some formulas from our previous paper [17]. In the /^-factorization ap- 
proach, the differential cross section for inclusive heavy quark production may be written 

as 

da(pp -^QQX) = —— ^ ql T , /i 2 ) $(x 2 , q 2T , /i 2 ) x 

ioir{XiX 2 s) , , 

x 2^ I M Isha(^ 5 -»■ QQ) dyi dy 2 dp 2T dq 1T dq 2T — — — , 

where g 2 T , /i 2 ) and $(x 2 , q 2T , /i 2 ) are the unintegrated gluon distributions in the proton, 

0i, 02 and 0q are the azimuthal angles of the initial gluons and final heavy quark respectively, 
Yl \M\snA{g*g* — > QQ) is the off-mass-shell matrix element. Symbol in Equ. (7) indicates 
the averaging over the initial and the summation over the final polarization states. The 
expression for the ^ |M|5 HA (g*g* — > QQ) was obtained in our previous paper [17]. 

Formulas for the differential cross sections in the standard parton model (SPM) may 
be obtained from Equ. (7) if we take the limit q\ 2T — > and average over the transverse 
momentum vector qi t2 r' 

dcr(pp^ QQX) = 1 —xiG(xi, fi 2 )x 2 G(x 2 , /i 2 ) x 

IQit(xiX 2 s) 2 ( v 

x \ M \lu(99 -> QQ) dj/i ^2 ^it 

where |M|p M (<7<? — > QQ) is the gluon-gluon fusion matrix element obtained in the standard 
parton model. 



2.3 Unintegrated gluon distributions 

Various parametrizations of the unintegrated gluon distribution used in our calculations 
are discussed below (see also [23]). 

As the first set, we use a BFKL-like parametrization (hereafter denoted as the JB 
parametrization) given in [24]. The method proposed in [24] lies upon a straightforward 
perturbative solution of the BFKL equation where collinear gluon density xG(x, /i 2 ) is used 
as the boundary condition. The unintegrated gluon distribution is calculated as a convolu- 
tion of collinear gluon distribution xG(x, /i 2 ) with universal weight factors: 

i 

$(x,g 2 ,/i 2 ) = J <p(ri,q^,fj, 2 )^G f/j drj, (9) 



X 



where 



as 

-,2 



Jo^asHl/^ln^/q^, if q\ < f 



V (V,q 2 T ^ 2 ) = { r] ^ } . < (10) 

^ J„ \2yfas ln(l/ V ) ln(^)j , if g 2 > »\ 

where Jo and Jo stand for Bessel functions of real and imaginary arguments respectively, 
and a s = 3as/n. In calculations we used the standard GRV set [25] for xG(x,fj, 2 ). The 
parameter as is connected with the Pomeron trajectory intercept: A = 4o,sln2 in the LO, 
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and A = 4as ln2 — Na 2 s in the NLO approximations, where N ~ 18 [26, 27]. The latter value 
of A could have dramatic consequences on the high energy phenomenology. However, some 
resummation procedures proposed in the last years lead to positive values: A ~ 0.2 — 0.3 [27, 
28]. The result A = 0.35 was obtained from the description of the px spectrum of D* mesons 
in the electroproduction at HERA [29]. We use this value of the parameter A in the present 
paper 4 . 

Another set of unintegrated gluon densities (the KMS parametrization) [30] is obtained 
from a unified BFKL and DGLAP description of F 2 data and includes the so called con- 
sistency constraint [31]. The consistency constraint introduces a large correction to the 
LO BFKL equation: about 70% of the full NLO corrections to the BFKL exponent A are 
effectively included in this constraint, as is explaned in [32]. 

Finally, the third unintegrated gluon function used here is the one which is obtained from 
conventional gluon density xG(x,fj, 2 ) by taking the /^-derivative [7, 11, 33]: 

»(*,,>) = ££M . (ii) 



=<?T 



Here we have used the expression for xG(x,[i 2 ) from the standard GRV set [25]. We point 
out that the parametrization (11), in contrast with JB and KMS parametrizations, takes into 
account the terms \n n (fi 2 / Aq CD ) and In™ (/i 2 / Aq CD ) ln n (l/:r) only. It is interesting to 
compare the numerical results obtained with the parametrizations JB and (11) because they 
are derived from the same collinear density but underwent evolution according to different 
(BFKL or DGLAP) equations. 

The integration limits in (7) and (8) are given by the kinematic conditions of the D0 and 
CDF experiments [1-5]. The calculation of the heavy quark production cross section in the 
fey-factorization approach have been done according to (7) for q\ T > Q 2 , and q\ T > Ql- For 
the regions q\ T < Qq and q\ T < Q 2 ,, we set q\ T = and q\ T = in the matrix elements of 
the hard subprocesses, take Yl \M\p M (gg — > QQ) instead of Yl \M\l UA (g*g* — > QQ) and use 
equation (8) of the usual parton model. The contributions from the assymmetric configura- 
tions (qf T < Qq, q\ T > Qq and q\ T > Qq, q\ T < Qq) are included in a similar way, where one 
of the gluons is described by the unintegrated distribution and the other one by the collinear 
density. The choice of the critical value of the parameter Ql — 1 GeV 2 is determined by 
the requirement that the value of as(^ 2 ) in the region q\ 2 r — Qo be small (where in fact 
a s (fi 2 ) < 0.26). 



3 Numerical results 

In this section we present the numerical results of our calculations and compare them 
with the D0 [4, 5], CDF [1-3] and UA1 [34] data. 

Besides the choice of the unitegrated gluon distribution, our theoretical results depend on 
the bottom quark mass, the factorization scale fj, 2 and the b quark fragmentation function. 
For example, a special choice of the 6-quark fragmentation function was used in paper [6] as 
a way to increase the 5-meson production cross section in observable region of transverse 
momentum. In the present paper we convert b quarks into 5-mesons using the usual Peterson 

4 We also used this value of A in the analysis of experimental data on the J/ip photo- and lcptoproduction at HERA [19, 20] 
and in the description [21, 22] of deep inelastic structure functions F£, F£ and in the small x region. 
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fragmentation function [35] with e = 0.006. Regarding the other parameters, we use m& = 
4.75 GeV and /i 2 = q\ 2T as in [8, 14] 5 . 




5 6 7 3 9 10 11 12 13 H 15 



PTmin- GeV 



Figure 2: The b quark transverse momentum distributions at Tevatron conditions presented in the form 
of integrated cross sections. The cuts applied: \yi\ < 1.5, 1 2/2 1 < 1.5, y/s — 630GeV (Fig. 2a), \yi\ < 1, 
y/s = 1800 GeV (Fig. 2b) and \yi\ < 1, \y 2 \ < 1, a/s = 1800 GeV (Fig. 2c). Curve 1 corresponds to the SPM 
calculations in the leading order approximation with GRV (LO) gluon density, curves 2, 3 and 4 correspond 
to the fcT-factorization results with JB, KMS and Equ. (11) unintegrated gluon distributions. Experimental 
data are from UAl [34] A, D0 [5] • and CDF [3] ■. 

The results of our calculations are shown in Fig. 2 — 9. Fig. 2 displays the b quark 
transverse momentum distribution at Tevatron conditions presented in the form of integrated 
cross sections. The following cuts were applied: \yi\ < 1.5, 1 2/2 1 < 1-5, = 630 GeV 
(Fig. 2a), \ yi \ < 1, y/s = 1800 GeV (Fig. 2b) and |yi| < 1, \y 2 \ < 1, y/a = 1800 GeV 
(Fig. 2c). Curve 1 corresponds to the SPM calculations at the leading order approximation 
with the GRV (LO) gluon density, curves 2, 3 and 4 correspond to the fc^-factorization 
results with the JB, KMS and the (11) unintegrated gluon distributions, respectively. One 
can see that the results obtained in the /cy-factorization approach agree very well with the 
D0 and CDF experimental data. The calculations based on the parametrization (11), which 
takes into account only the terms ln n (/i 2 /AQ CD ) and \n n (ft 2 / Aq CD ) ln n (l/x) predict 
the cross section which is lower than the data by a factor of about 2 (Fig. 2a and Fig. 2b). 
We would like to note the difference in the shapes between the curves obtained using the 

5 We also used this choice of fi 2 earlier [19, 20] for the description of J/ip photo- and leptoproduction processes at HERA 
within the fcy-factorization approach and the colour singlet model. 
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/c-r-factorization approach and the standard parton model. This difference shows the pr 
broadening effect which is usual for the /cy-factorization approach. At the same time, the 
shape of the curves 4 and 1 is practically identical. 




Figure 3: The ratio of a(b) at s/s = 630 GcV to s/s = 1800 GeV as a function of the minimum b quark 
transverse momentum Px m in- Notation of the curves 1 and 3 is the same as in Fig. 2. Experimental data 
are from CDF[1] ■ 




5 10 15 20 25 
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Figure 4: Theoretical predictions for the B meson px spectrum at y/s = 1800 GeV compared to the CDF 
data. The cuts applied: \y B \ < 1. Notation of the curves 1 — 4 is the same as in Fig. 2. Experimental data 
are from CDF [2] ■. 

It is notable that, in the collinear approximation, the sum of the LO and NLO pQCD con- 
tributions still underestimates the b quark production rate by a factor of 2 [3, 5]. The results 
obtained in the fc^-factorization approach in Ref. [16] lie above the sole LO contribution, but 
below the sum of LO and NLO contributions, in an apparent disagreement with our present 
results. The roots of this discrepancy are connected with the parameter settings accepted in 
Ref. [16], that is, the large values of the quark mass m& = 5 GeV and the remormalization 
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scale /i 2 = ml 2T in the running coupling constant. 




Figure 5: The cross section for muons from B meson decay as a function of rapidity compared to the D0 
data. The cuts applied: p^, > 5 GeV (Fig. 5a) and p£ > 8 GeV (Fig. 5b) . Notation of the curves 1 — 4 is 
the same as in Fig. 2. Experimental data are from D0 [4] •. 

Fig. 3 shows the ratio of the cross sections measured at different beam energies, a(b) 
at yfs = 630 GeV to y/s = 1800 GeV, as a function of the minimum b quark transverse 
momentum p^ min . Notation of the curves 1, 3 is the same as in Fig. 2. One can see that the 
experimental data collected by the D0 collaboration agree with the LO pQCD calculations 
as well as with the fc^-factorization ones. This result is not surprising because, when the 
ratio of the cross sections is considered, many factors affecting the absolute normaliszation, 
as well as many theoretical and experimental uncertainties partially or completely cancel 
out [1]. 




|p T '|, GeV 



Figure 6: Predictions on the leading muon px spectrum in the bb production events compared to the D0 
data. The cuts applied to both muons: 4 < p% < 25 GeV, \^\ < 0.8 and 6 < m w < 35 GeV. Notation of 
the histograms 1 — 4 is the same as in Fig. 2. Experimental data are from D0 [5] •. 
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Fig. 4 shows the prediction for the B meson px spectrum at = 1800 GeV compared 
to the CDF data [3] within the experimental cuts \y B \ < 1. Notation of the curves 1 — 4 
is the same as in Fig. 2. Here we find good agreement between the results obtained in the 
/cy-factorization approach and experimental data. One can see that the pt broadening effect 
mentioned earlier appears to be not as much clear for the B meson production as it was 
for the b quark production. This is because of compensatory effects in the fragmentation 
process. We note again that the ordinary NLO pQCD calculations underestimate the B 
meson production by a factor of 3 [2]. 




Figure 7: Double differential cross section for muons from B meson decay as a function of pi^ compared to 
the D0 data [4]. The cuts applied: 2.4 < \y^\ < 3.2. Notation of the curves 1 — 4 is the same as in Fig. 2. 
Experimental data are from D0 [4] •. 

The recent D0 experimental data refer also to muons originating from semileptonic decays 
of 5-mesons. To produce muons from 5-mesons in theoretical calculations, we simulate 
their semileptonic decay according to the standard electroweak theory. In Fig. 5 we show 
the rapidity distribution da/dly^l of the decay muons for both pt^ > 5 GeV (Fig. 5a) and 
Pt > 8 GeV (Fig. 5b). Notation of the curves 1 — 4 is the same as in Fig. 2. We find that 
the /cr-factorization results agree very well with the D0 experimental data. In the central 
rapidity range \y^\ < 1, the results based on the parametrization Equ. (11) lie below the 
ones based on JB and KMS parametrizations by a factor of about 1.5. The ordinary NLO 
pQCD calculations lie below the data by a factor of about 4 [4]. 

Fig. 6 shows the leading muon p? spectrum in the bb production events compared to the 
D0 data, where the leading muon is defined as the muon with the greatest pj. value. The 
cuts applied to both muons are 4 < p^, < 25 GeV, < 0.8 and 6 < < 35 GeV. 
Notation of the histograms 1 — 4 is the same as in Fig. 2. One can see that the histograms 
2 and 3 lie even a bit higher than the experimental data. 

The double differential cross sections da/dp^dy^ 1 in the forward rapidity region 2.4 < 
\y^\ < 3.2 are also well described by the /^-factorization approach (Fig. 7). Notation of the 
curves 1 — 4 is the same as in Fig. 2. It is interesting to note that the shape of all the 
curves is practically the same. The NLO pQCD calculations underestimate the D0 data by 
a factor of 4 [4]. 

We point out that the investigations of bb correlations such as the azimuthal opening 
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Figure 8: Azimuthal muon-muon correlations at Tevatron conditions. The cuts applied to both muons: 
4 < < 25 GeV, < 0.8 and 6 < m w < 35GeV. Notation of the histograms 1 — 4 is the same as in 
Fig. 2. Experimental data are from D0 [4] •. 

angle between the b and b quarks (or between their decay muons) allow additional details 
of the b quark production to be tested since these quantities are sensitive to the relative 
contributions of the different production mechanisms [8, 12-14, 16]. In the naive gluon- 
gluon fusion mechanism, the distribution over the azimuthal angle difference A(p bb must be 
simply a delta function 5(A(f> bb — n). Taking into account the non- vanishing initial gluon 
transverse momenta qy? and q2T leads to the violation of this back-to-back quark production 
kinematics in the /^-factorization approach. 

The differential bb cross section da/dAcj)^ is shown in Fig. 8. The cuts applied to both 
muons are 4 < pt^ < 25 GeV, |r/ M | < 0.8 and 6 < m w < 35 GeV. Notation of the histograms 
1 — 4 is the same as in Fig. 2. One can see that good agreement between JB and KMS 
predictions and the experimental data is observed. The shape of the histogram 4 strongly 
differs from that of the histograms 2 and 3. In the small Acf)^ ~ region, the parametrization 
(11) underestimates the D0 experimental data. This fact indicates the importance of the 
large ln n (l/a;) contributions. One can see that the properties of different unitegrated 
gluon distributions manifest themselves in the bb or muon-muon azimuthal correlations. As 
expected, the sole LO pQCD contribution predicts a peak at A0 W ~ n. 

In addition, we present our predictions for muon-muon (Fig. 9a) and muon-jet (Fig. 9b) 
cross sections at the Tevatron and CERN LHC conditions. In the latter case, we take 
the kinematic requirements of the detector ATLAS [36] as a representative example. This 
conditions imply the presence of a decay muon with > 6 GeV and < 2.5. Curves 1 — 4 
are the same as in Fig. 2. We point out that the predicted cross sections are rather uncertain 
and approximate because we neglected the saturation effects in the gluon distributions [8, 
12, 13]. However, the treshold value of x (where these effects will come into play) is still 
unknown. 
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Figure 9: Theoretical predictions for muon-muon (Fig. 9a) and muon-jet (Fig. 9b) cross sections as functions 
of y/s at the Fermilab Tevatron and CERN LHC conditions. The cuts applied: pt^ > 6 GeV, |y M | < 2.5. 
Notation of the curves 1 — 4 is the same as in Fig. 2. 

4 Conclusions 

In this paper we have considered heavy quark production in pp collisions at Tevatron in 
the framework of the /cr-factorization approach. We investigated the dependence of the b 
quark, B meson and the decay muon cross sections on different forms of the unintegrated 
gluon distribution. The analysis covered the azimuthal correlations between the b and b 
quarks and their decay muons. We compared the theoretical results with recent experi- 
mental data collected by the D0 and CDF collaborations at Tevatron. We found that the 
/^-factorization results agree well with the experimental data when we use JB or KMS un- 
integrated gluon distributions and set the quark mass = 4.75 GeV, the factorization scale 
fi 2 = and Aq C d = 250 MeV. The properties of different unitegrated gluon distributions 
manifest themselves in the bb or muon-muon azimuthal correlations. From the analysis of 
these correlations we can conclude that JB and KMS unintegrated gluon distributions are 
more preferable than the parametrization (11). Finally, we present our predictions for the 
muon-muon and muon-jet cross sections at Tevatron and CERN LHC conditions. 



5 Acknowledgments 

The study was supported in part by RFBR grant N° 02-02-17513. A.V.L. also was 
supported by INTAS grant YSF 2002 N° 399. 

References 

[1] D. Acosta et al. (CDF Collab.), hep-ex/02060191 

[2] D. Acosta et al. (CDF Collab.), Phys. Rev. D 65, 052005 (2002). 



11 



F. Abe et al. (CDF Collab.), Phys. Rev. D 55, 2546 (1997). 
B. Abbott et al. (D0 Collab.), hep-ex/9907029| 

B. Abbott et al. (D0 Collab.), Phys. Lett. B 487, 264 (2000); |hep^"ex/9905024| 
M. Cacciari and P. Nason, |hep-ph/0204025| 

L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100, 1 (1983). 

E.M. Levin, M.G. Ryskin, Yu.M. Shabelsky and A.G. Shuvaev, Yad. Fiz. 53, 1059 
(1991). 

S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B 366, 135 (1991). 
J.C. Collins and R.K. Ellis, Nucl. Phys. B 360, 3 (1991). 

E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 44, 443 (1976); 45, 199 
(1977); 

Yu.Yu. Balitsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 
M.G. Ryskin and Yu.M. Shabelski, Z. Phys. C 69, 269 (1996). 

M.G. Ryskin, Yu.M. S habelski and A.G . Shuvaev, Phys. Atom. Nucl. 64, 1995-2005 
(2001); hep-ph/990757; |hep-ph/0007238| 

P. Hagler, R. Kirschner, A. Schafer et a/., Phys. Rev. D 62, 071502 (2000). 

H. Jung, Phys. Rev. D 65, 034015 (2002); |hep^ph/0110034| 

S.P Baranov and M. Smizanska, Phys. Rev. D 62, 014012 (2000). 



in 



A.V. Lipatov, V.A. Saleev and N.P Zotov, Yad. Fiz. 66, 1 (2003); hep -ph/0112114 
Proceedings of the QFTHEP'2001, Moscow, Russia, 2002, p. 238. 

V.A. Saleev and N.P. Zotov, Mod. Phys. Lett. A 9, 151 (1994); A 11, 25 (1996); 

A.V. Lipatov and N.P. Zotov, Mod. Phys. Lett. A 15, 695 (2000); 

A.V. Lipatov, V.A. Saleev and N.P. Zotov, Mod. Phys. Lett. A 15, 1727 (2000). 

A.V. Lipatov and N.P. Zotov, to be published in Yad. Fiz. (2003); |hlp~ph/0208237| 

A.V. Lipatov and N.P. Zotov, to be published in Eur. Phys. J. C (2003); |hep^ph/0210310| 

A.V. Kotikov, A.V. Lipatov, G. Parente and N.P. Zotov, Eur. Phys. J. C 26, 51 (2002). 

A. V. Kotikov, A .V. Lipatov and N.P. Zotov, to be published in Eur. Phys. J. C (2003), 
|hep-ph/0207226l 

B. Andersson et al. (Small x Collab.), Eur. Phys. J. C 25, 77 (2002). 
J. Bliimlein, DESY 95-121. 

M. Gluck, E. Reya and A. Vogt, Z. Phys. C 67, 433 (1995). 

V.S. Fadin and L.N. Lipatov, Phys. Lett. B 429, 127 (1998); 
M. Ciafaloni and G. Camici, Phys. Lett. B 430, 349 (1998). 



12 



[27] G. Salam, JHEP 9807:019 (1998); |hep^p h/9806482, 

[28] S.J. Brodsky, V.S. Fadin, V.T. Kim, L.N. Lipatov and G.B. Pivovarov, JETP Lett. 70, 
155 (1999). 

[29] S.P Baranov and N.P Zotov, Phys. Lett. B 458, 389 (1999); B 491, 111 (2000). 

[30] J. Kwiecinski, A.D. Martin and A.M. Stasto, Phys. Rev. D 56, 3991 (1997). 

[31] J. Kwiecinski, A.D. Martin and A.M. Sutton, Phys. Rev. D 52, 1445 (1995); Z. Phys. 
C 71, 585 (1996). 

[32] J. Kwiecinski, A.D. Martin and J.J. Outhwaite, Eur. Phys. J. C 9, 611 (2001). 
[33] N.N. Nikolaev and B.G. Zakharov, Phys. Lett. B 333, 250 (1994). 
[34] C. Albajar et al, Z. Phys. C 61, 41 (1994). 

[35] K.A. Peterson, D. Schlatter, I. Schmitt and P.M. Zerwas, Phys. Rev. D 27, 105 (1983). 
[36] W. Armstrong et al. (ATLAS Collab.), CERN/LHCC/94-43. 



13 



